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Abstract 

O ! 

Understanding the transport properties of a porous medium from a knowledge of its microstruc- 
ture is a problem of great interest in the physical, chemical and biological sciences. Using a first- 

<n : 

\Q , passage time method, we compute the mean survival time r of a Brownian particle among perfectly 

(N ■ 

absorbing traps for a wide class of triply-periodic porous media, including minimal surfaces. We 
find that the porous medium with an interface that is the Schwartz P minimal surface maximizes 

OS : 

the mean survival time among this class. This adds to the growing evidence of the multifunctional 

. ^ ' optimality of this bicontinuous porous medium. We conjecture that the mean survival time (like 

>< ■ 

the fluid permeability) is maximized for triply periodic porous media with a simply connected 

a ; 

pore space at porosity <fi = 1/2 by the structure that globally optimizes the specific surface. We 
also compute pore-size statistics of the model microstructures in order to ascertain the validity of 
a "universal curve" for the mean survival time for these porous media. This represents the first 
nontrivial statistical characterization of triply periodic minimal surfaces. 
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I. INTRODUCTION 



Fluid-saturated porous media are ubiquitous in nature (e.g., geological and biological 
media) and in synthetic situations (e.g., filters, cements and foams) Understanding 
the transport properties of a porous medium from a knowledge of its microstructure is a 
subject that spans many fields , |6j. In particular, physical phenomena involving 

simultaneous diffusion and reaction in porous media abound in the physical and biological 
sciences , |7| . Considerable attention has been devoted to instances in which diffusion 
occurs in the pore region of the porous medium with a "trap" region whose interface can 
absorb the diffusing species via a surface reaction. Examples are found in widely different 
processes, such as heterogeneous catalysis, fluorescence quenching, cell metabolism, ligand 
binding in proteins, migration of atoms and defects in solids, and crystal growth A key 
parameter in such processes is the mean survival time r, which gives the average lifetime of 
the diffusing species before it gets trapped. 

It is noteworthy that while there has been a significant amount of progress made on 
the determination of the structures that optimize a variety of transport and mechanical 
properties of porous media [l|, s|, Q], there have been no studies that have attempted to 
find the optimal isotropic structures for the mean survival time r [h]]. It has been recently 
demonstrated that triply periodic two-phase bicontinuous composites with interfaces that 
are the Schwartz P and D minimal surfaces are not only geometrically extremal, but are also 
extremal for simultaneous transport of heat and electricity [ll, I2I, 13 1 . A minimal surface 
is one that is locally area minimizing i.e., every point has zero mean curvature. Triply 
periodic minimal surfaces are minimal surfaces that are periodic in all three coordinate 
directions. An important subclass of triply-periodic minimal surfaces are those that partition 
space into two disjoint but intertwining regions that are simultaneously continuous (i.e., 



bicontinuous) 14], |l5j. Examples of such surfaces include the Schwartz primitive (P), the 



Schwartz diamond (D), and the Schoen gyroid (G) surfaces [see Fig. [T]; each disjoint region 
has volume fraction equal to 1/2. Such triply periodic minimal surfaces arise in a variety 

is ll7L block copolymers [3] as 
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22]. 



of systems, including nanocomposites [16(, micellar material 
well as lipid-water systems and certain cell membranes []i 

The multifunctionality of such two-phase systems has been further established by show- 
ing that they are also extremal when a competition is set up between the effective bulk 
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FIG. 1: (Color online) Unit cells of three different minimal surfaces, (a) Schwartz P surface, (b) 
Schwartz D surface, (c) Schoen G surface. Image is adapted from Ref. [24]. 



modulus and electrical (or thermal) conductivity of the bicontinuous composite [23j. Jung 



and Torquato 



241 ] also explored the multifunctionality of the three minimal surfaces shown 



in Fig. [I] with respect to Stokes (slow viscous) flow. The simulations conducted by Jung and 
Torquato revealed that the Schwartz P porous medium has the largest fluid permeability k 
among a class of structures studied with porosity = 1/2. Further, the fluid permeability 
was found to be inversely proportional to the specific surface s (interface area per unit vol- 
ume). This led the authors to conjecture that the maximal fluid permeability for a triply 
periodic porous medium with a simply connected pore space [25j at a porosity 0=1/2 is 



achieved by the structure that globally minimizes the specific surface 

In this manuscript, we explore whether the mean survival time r of a Brownian particle 
among perfectly absorbing traps is maximized by the Schwartz P structure among a wide 
class of triply-periodic porous media based on the class of six models studied in Ref. 24]. 
Further, we go beyond considering just a representative medium from each class of mod- 
els and explore several of the models over a range of parameter values, as will be further 
explained below. This inquiry is also motivated by a cross-property bound that rigorously 
links the fluid permeability k to the mean survival time r 261 ]: 



k < r. 



(1) 



It should be noted that the inequality of ([I]) becomes an equality for transport in parallel 
cylindrical tube bundles of arbitrary cross-sectional geometry 26l |27|. It is also a relatively 
tight bound (when appropriately scaled) for transport around distributions of inclusions 
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261 ] . It is clear that the permeability k will be maximized for a simply connected pore space 



(e.g., the presence of dead ends are undesirable because they would not contribute to the 



flow). The inequality ([T]) and the results of Ref. [24| suggest that r is maximized by the 
same simply connected microstructure that maximizes k for <fi = 1/2. 

We will also test whether our calculations for r collapse on to a "universal" scaling relation 
for the mean survival time 281 ] . This requires us to compute the pore-size density functions 
for the triply periodic surfaces considered in this paper. 

In Sec. II, we define terminologies and give a precise statement of the problem. In 
Sec. Ill, we describe the first-passage time technique that we utilize to compute the mean 
survival time for a wide class of triply periodic porous media that are generally bicontinuous. 
Section IV reports our finding for r. In Sec. V, we discuss a universal scaling relation for r, 
report pore-size statistics for the triply periodic porous media studied here, and ascertain 
the applicability of the universal curve for these structures. Finally, in Sec. VI, we make 
concluding remarks and discuss the ramifications of our results. 

II. DEFINITIONS AND PROBLEM STATEMENT 

The mean survival time r arises in steady-state diffusion of reactants in a trap-free pore 
region Vi with diffusion coefficient D among static traps with a unit rate of production of the 
reactants per unit pore volume [lj. When the reactants come in contact with the pore-trap 
interface <9V, they get absorbed. Using homogenization theory [l], it has been shown that 
r, the average time traveled before a diffusing particle gets trapped, is given by: 

r=Q, CO 



where u(r) is the scaled concentration field of reactants, which satisfies the diffusion equation 

Au = -1 in Vi (3) 
u = on dV, (4) 

and A is the Laplacian operator. 



Following Ref. 24] , we compute r for a wide class triply periodic structures: the Schwartz 
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P, Schwartz D and Schoen G minimal surfaces [Fig. [T], a cubic pore-square channel model 
[Fig. [2(a)], and a spherical pore-circular channel model [Fig. [21(b)], and an array of spherical 
traps arranged on a simple cubic lattice Q] (which of course is not bicontinuous) . The 
pore-channel models each contain two parameters, a and b, that enable one to control the 
relative size of the pores (either spherical or cubic in shape) to the size of the channels (either 
circular or square cross-sections) 24J. The parameter a determines channel width/diameter, 
whereas b determines the width/diameter of the pore. This class of models is bicontinuous 
provided that a > 0. In the current investigation, we compute r for the pore-channel models 
using a wider range of values for a and b than used in Ref. {24], while keeping fixed at 
1/2. 
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FIG. 2: (Color online) Unit cells of two different pore-channel models, (a) Cubic pore-square 
channel and its cross-section, (b) Spherical pore-circular channel and its cross section. Image is 
adapted from Ref. [24]. 
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III. FIRST-PASSAGE TIME METHOD FOR COMPUTING MEAN SURVIVAL 
TIME 



The mean survival time r for the aforementioned triply periodic porous media could 
be calculated using standard random-walk techniques that simulate the detailed zig-zag 
trajectory of a Brownian particle {^J. We instead compute r for these microstructures 



using an efficient first-passage time algorithm 31] • The fundamental periodic cell is taken 
to be a cube of side length unity. 

The first-passage time algorithm is implemented by applying the following set of rules: 

1. Introduce a Brownian particle into a random position in the trap-free phase. 

2. While the walker is sufficiently far from the two-phase interface, construct the largest 
sphere of radius Ri centered at the Brownian particle that just touches the two-phase 
interface. 

3. In one step, the walker jumps to a random point on the surface of the sphere [Fig. [3], 
with an average hitting time of: 

i=R 2 J{QD). (5) 

4. Repeat steps (2)-(3) until the walker is within some prescribed small distance 5 (taken 
to be 10~ 8 in our simulations) from the two-phase interface. 

5. Repeat steps (l)-(4) for many random walkers and calculate the mean survival time 
from the following equation: 

where angular brackets denote an ensemble average. 

The mean survival time r of the pore-channel models are determined using this algorithm. 
However, since we only have discrete (not continuous) representations of the minimal sur- 



faces [24j, the discrete analog of the algorithm is utilized in these cases |32j. In the discrete 
case, first-passage cubes of length 2Ri are utilized instead of first-passage spheres, and the 
average time it takes to move to the surface of a first-passage cube is t « 0.225i?f /D. 



FIG. 3: (Color online) Two-dimensional depiction of continuous first-passage time method applied 
to the spherical pore-circular channel model. The Brownian particle walks, by jumping to the 
surface of the largest possible first-passage sphere at each step, until it gets trapped at the two- 
phase interface (maroon boundary layer, which appears as dark grey if the image is being viewed 
in black and white). 

IV. RESULTS FOR THE MEAN SURVIVAL TIMES 

The value of r for each of the aforementioned models is given in Table HJ along with the 
fluid permeability and specific surface. All fluid permeability measurements (if available), 
as well as the specific surface measurements for the minimal surfaces, and the pore-channel 

n 

models with 6 = are taken from Ref . [24( . We see that r always obeys the rigorous bound 
specified by (1). We also see that a simply connected pore phase is a crucial topological 
feature required to achieve large mean survival times at a porosity <fi = 1/2. However, a 
simply connected pore phase is not a sufficient condition, as evidenced by the relatively 
small mean survival times associated with the Schwartz D and Schoen G minimal surfaces. 
These two minimal surfaces have rather large specific surfaces and hence serve as efficient 
traps for the diffusing Brownian particles. Indeed, we see that the survival times are inversely 
proportional to the corresponding specific surfaces for all of the bicontinuous structures. The 
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pore spaces of structures with large r are expected to be simply connected and therefore it 
would not be unreasonable for the survival time to be inversely proportional to the specific 
surface in these instances. As hypothesized, we find that among the structures examined, 
the Schwartz P porous medium maximizes the mean survival time. To get some idea of 
the fluctuations about the average r values, we computed the associated variance a\ for 
trapping for the three minimal surfaces and found that cr^ = 0.00051,0.00013 and 0.00005 
for the Schwartz P, Schoen G and Schwartz D surfaces, respectively. 



TABLE I: Mean survival time r, fluid permeability k, and specific surface s of triply periodic 
structures. All quantities are made dimensionless using the side length of the unit cell. 



Structure 


r 


k 


s 


Schwartz P 


0.0173950 


0.0034765 


2.3705 


Schoen G 


0.0093266 


0.0022889 


3.1284 


Schwartz D 


0.0060414 


0.0014397 


3.9011 


Cubic-pore channel 








(a = 0.25; 6 = 0) 


0.0139289 


0.0030744 


3.0000 


(a = 0.1324; 6 = 0.25) 


0.0123813 


0.0005310 


3.8360 


(a = 0.0781; 6 = 0.3125) 


0.0122259 


0.0000948 


3.9254 


(a = 0; 6 = 0.3969) 


0.0127852 




3.7804 


Sphere-pore channel 








(a = 0.2836; 6 = 0) 


0.0167934 


0.0034596 


2.63990 


(a = 0.1480; 6 = 0.2794) 


0.0164995 




2.73649 


(a = 0.0908; 6 = 0.3633) 


0.0161733 




2.93057 


(a = 0; 6 = 0.4924) 


0.0161093 




3.04681 


Spherical trap (0 = 0.5) 


0.0139640 


0.0030591 


3.0780 



Given that both the mean survival time r and fluid permeability k are made dimensionless 
with the side length of the unit cell, one might ask why these dimensionless bulk properties 
are two to three orders of magnitude smaller than unity? This is a well-known behavior 
for such transport properties of porous media, optimal or not . For example, the fluid 
permeability (average of the velocity field) can be regarded to be the "effective pore channel 
area of the dynamically connected" part of the pore space. For non-simply connected pore 
spaces, there generally will be regions that contain fluid but do not actively contribute 
to the flow (not dynamically connected). Moreover, because of the no-slip condition, the 
velocity only becomes significantly large sufficiently away from the pore-solid interface for 
general porous media. Thus, these two effects conspire to make the effective area of a pore 
channel, i.e., the fluid permeability k, considerably less than the geometric pore channel 



S 



sizes. (This effect is even true for the simple case of flow in a tube; see Ref. |27j). Similar 
arguments apply to the mean survival time. The perfectly absorbing boundary condition 
at the pore-solid interface (i.e., zero concentration field) means that the concentration field 
becomes appreciably large sufficiently away from the pore-solid interface, which results in 
a mean survival time (average of the concentration field) that is generally several orders of 
magnitude smaller than dictated by the largest pore dimensions. 

V. PORE-SIZE FUNCTIONS AND UNIVERSAL SCALING RELATION FOR 
SURVIVAL TIME 

A. Pore-Size Function 

Porous media whose interfaces are triply periodic minimal surfaces apparently have re- 
markable macroscopic properties. Nonetheless, these structures have yet to be statistically 
characterized using nontrivial descriptors. Here we present pore-size functions for these 
structures as well as the triply periodic circular channels. This is motivated by the fact 
that the pore-size density function P(S) arises in rigorous lower bounds on r [l| as well as a 



universal curve for r [281] . described below. The quantity P(5)d5 gives the probability that 
a randomly chosen point in the pore region lies at a distance between 5 and 5 + d5 from the 
nearest point on the interface jj. 

We employ the following algorithm to determine the pore-size function P(S) [331 ]: 

1. Choose a random location in the pore phase. 

2. Find the radius of the largest sphere centered at the above point that just touches the 
two-phase interface. 

3. Repeat steps 1-2 for many random locations and create a list of radii. 

4. After sampling sufficiently, bin the sphere radii. Divide the number of radii in each 
bin by the total number of radii to determine P(5). 

The interpolated functional form of P(8) is given in Eq. (j7j), flS}, flH]) and (fTOl) for the Schwartz 
P, Schwartz D, Schoen G minimal surfaces, and the circular channel model (spherical pore 
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circular-channel model with 6 = 0), respectively [see Fig. Hj: 

, -1495 5 + 1665 4 - 61.25 3 + 8.25 2 - 0.6965 + 0.164, if 5 < 0.419 
P(S) ={ (7) 

0, otherwise 



, 22385 s - 5235 4 - 44.15 3 + 12.85 2 - 1.555 + 0.275, if 5 < 0.202 
P(5) = { (8) 

0, otherwise 



, -13885 5 + 8975 4 - 2255 3 + 23.45 2 - 1.475 + 0.223, if 5 < 0.234 
P(5) = { (9) 

0, otherwise 



. -1625 5 + 1735 4 - 59.95 3 + 7.575 2 - 0.7645 + 0.176, if 5 < 0.234 
P(5) = { (10) 

0, otherwise. 

For the spherical pore (spherical pore circular-channel model with a = 0), the pore size 
density function is known exactly: 



, if5<a , 

P(5) = { a 3 (11) 

0, otherwise. 



B. Universal Curve 

Based on rigorous lower bounds on the survival time [l| , the following "universal" scaling 
relation for r has been found to apply to a wide class of microstructures and range of 
porosities [28] : 

- = -x + - X \ (12) 
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FIG. 4: (Color online) Pore size density function data generated from the presented algorithm 
for the three minimal surfaces and the circular-channel model, along with the best-fit fifth degree 
polynomial. 

where r = 302/(-D0s 2 ), x = (5) 2 /(toD), and (5) is the mean pore size, denned by the first 
moment of the pore-size density function P(5) 



(5) = I 5P(8) dS. 
'o 



(13) 



To explore the robustness of the universal curve, we check if the reported mean survival 
time data for a subset of the triply periodic surfaces falls on the curve. In order to compare 
the reported r results to Eq. (112j) . we determine (5) from the pore-size density function 
P(S) results given above. Figure reveals that the scaled survival times of the minimal 
surfaces fall relatively close to the universal curve. The corresponding values for the other 
bicontinuous structures considered here are close to the minimal-surface values and hence 
are not shown in the figure. 
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FIG. 5: (Color online) The dimensionless mean survival time t/tq versus the scaled mean pore size 
squared, (5) 2 /(tqD), for various two-phase media. The solid curve is the universal scaling relation 
CEH) H- 



VI. CONCLUSIONS AND DISCUSSION 



To conclude, we have shown that a simply connected pore phase is a crucial topological 
feature required to achieve large mean survival times at a porosity = 1/2. However, 
a simply connected pore phase is not a sufficient condition feature; one must also have 
interfaces with small specific surface s. Indeed, we found that the survival times are inversely 
proportional to the corresponding specific surfaces for all of the triply periodic structures 
considered here. Nonetheless, the pore spaces of structures with large survival times are 
expected to be simply connected and therefore it would not be unreasonable for r to be 
inversely proportional to the specific surface in these instances. We have demonstrated, as 
hypothesized, that the Schwartz P porous medium maximizes the mean survival time among 
the set of twelve triply-periodic structures considered here. This lends further evidence to 
the multifunctional optimality of the Schwartz P minimal surface, making this structure 
of great practical value to guide the design of new materials with a host of desirable bulk 
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properties. Moreover, scaled survival times of the minimal surfaces fall relatively close to 
the universal curve ffl2l) found in Ref. [28 ] . 

Based on our findings, we conjecture that the mean survival time (like the fluid perme- 
ability) is maximized for triply-periodic porous media with a simply connected pore space 
at porosity = 1/2 by the structure the globally optimizes the specific surface. The verifi- 
cation of this conjecture remains an outstanding open question. This extremal problem falls 
in the general class of isoperimetric problems, which are notoriously difficult to solve. A pro- 
totypical isoperimetric example is Kelvin's problem: the determination of the space-filling 
arrangement of closed cells of equal volume that minimizes the surface area. Although it is 
believed that the Weaire-Phelan structure is an excellent solution to Kelvin's problem, 
there is no proof that it is a globally optimal one. Our conjecture is also likely a difficult 
one to prove. 

In this regard, it is noteworthy that an original goal of Ref. [l^| was to show that the 
triply-periodic surface with minimal specific surface s at porosity <p = 1/2 is the Schwartz P 
surface. While numerical simulations provided empirical evidence supporting this proposi- 



tion, the authors of Ref. [15] could not prove it rigorously. However, they were able to show 
that the Schwartz P, Schwartz D, and Schoen G minimal surfaces are local minima of the 
specific surface area s at fixed volume fraction = 1/2. Thus, the question of the global 
optimality of the Schwartz P surface (i.e., minimal total interface surface area or specific 
surface s), is an open question for future investigation. 
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